A novel diagrammatic technique for the single-site Anderson model 
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A diagrammatic theory around the atomic limit is proposed for the single-impurity Anderson 
model in which the strongly correlated impurity electrons hybridize with free (uncorrelated) con- 
duction electrons. Using this diagrammatic approach, we prove the existence of a linked cluster 
theorem for the vacuum diagrams and derive Dyson type of equations for the localized and conduc- 
tion electrons and corresponding equations for the mixed propagators. The system of equations can 
be closed by summing the infinite series of ladder diagrams containing irreducible Green's functions. 
The result allows to discuss resonances associated with the quantum transitions at the impurity site. 
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The study of strongly-correlated electron systems has 
become one of the most active fields of condensed matter 
physics during the last decade. The properties of these 
systems cannot be described by Fermi liquid theory. One 
of the important models of strongly correlated electrons 
is the single-site or impurity model introduced by Ander- 
son in 1961 [H, which has been intensively discussed in 
many papers. It is a model for a system of free conduc- 
tion electrons that interact with the system of a local spin 
of an electron in the d or f shells of an impurity atom. 
The impurity electrons are strongly correlated because of 
the strong on-site Coulomb repulsion and interact with 
the conduction electrons via exchange and hybridization. 
Most previous work is based on the method of equations 
of motion for retarded and advanced Green's functions 
and a truncating procedure as proposed, for example, by 
Bogoliubov and Tjablikov Q. A first attempt to develop 
a diagrammatic theory for this problem was undertaken 
by Barabanov in With introduction of the dynam- 
ical mean field theory, the interest in the Anderson im- 
purity model has considerably increased because infinite 
dimensional lattice models like the Hubbard model can 
be mapped onto effective impurity models and a self- 
consistency condition [1, (5[ . 

The Hamiltonian of the model taking only s-like elec- 
trons into account is written as 

H = H + Hi n t , Ho = Hq + H Q , 
H< = £e(k)c7+C k(7 , H£ = e f J2tiU + Unfn{, 

ka a 

H ™ = -^E(^/ ff + Ck CT + K CT c+ T / ff ), (1) 



where n s „ = f+f a , C k<7 (C£ CT ) and /<,(/+) arc the anni- 
hilation (creation) operators of conduction and impurity 
electrons with spin a, respectively; e(k) is the kinetic 
energy of the conduction band states (k, a); e/ is the lo- 



cal energy of /-electrons and N is the number of lattice 
sites. Hi n t is the hybridization interaction between con- 
duction and localized electrons. Summation over k will 
be changed to an integral over the energy e(k) using the 
density of states po(e) of conduction electrons. 

The term in the Hamiltonian involving U describes the 
on-site Coulomb interaction between two impurity elec- 
trons. This term is far too large to be treated by per- 
turbation theory. Therefore, we include it in the non- 
interacting Hamiltonian Hq. Since this term invalidates 
Wick's theorem for local electrons, we first of all have to 
formulate the generalized Wick's theorem (GWT) for lo- 
cal electrons, preserving the ordinary Wick theorem for 
conduction electrons. Our GWT can be considered as 
the prescription which allows to determine the irreducible 
Green's functions or Kubo cumulants. A similar prescrip- 
tion has been used when discussing the properties of the 
single-band Hubbard model @, 0, S| ■ 

In interaction representation, the renormalized (Mat- 
subara) Green's functions of conduction and impurity 
electrons have the form: 

G(k,a,r\ k>,V) =-{TC kff (T)C k v(r')^));, (2) 
g(a,T | a',r') = - <T/ ct (t)/ ct , {r')U(fi)) C Q . (3) 

Here r and r' stand for the imaginary time with < r < 
f3 (/? is the inverse temperature) and T is the chronologi- 
cal time ordering operator. The evolution operator U(/3) 
is determined by the hybridization interaction Hi nt . The 
statistical averaging in Eqs. (2) and (3) is carried out 
with respect to the zero-order density matrix of the con- 
duction and impurity electrons. 

The thermodynamic perturbation theory with respect 
to Hi n t requires an appropriate generalization scheme for 
the evaluation of statistical averages of the T-products of 
localized /-electron operators. 

In zero-order approximation (negelecting the hy- 
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bridization term), the corresponding Green's functions 
have the form (u = u> n — (2n + 1)tv/(3): 



G° aa , (k, k' | iu) = <W <5 kk , / {tw - e(k)} , 
1 - n w , n w 



g°a> = 5 acr ,g°(iu) = 5 a 



where (a = —a) and 



\<j(iuj) X w (iuj) 



(4) 



\ a {iuj) = iu} + E — E a , Xa{iuo) = iuj + E w - E 2 , 
Zo 

-f)E 2 



e -0E o + e -0E„ +e -(}Er + e -0E 2 



n w = {e-^ + e-^}/Z . 



In the case of /-electrons, we can formulate the identity 
corresponding to the GWT: 

(T/l/27 3 74> = ( T /l/4) ( T /2/3) 

- <T/i7 3 ) (Tf 2 7 4 ) + (T/i/ 2 / 3 7 4 )o r . ( 5 ) 



g 2 °(l,2|3,4) = 5 °(1|%°(2|3) 

- 5 °(l|3).9 (2|4)+.g(° )lr (l,2|3,4), (6) 

where n stands for (er n , r n ). The generalization for more 
complicated averages of type 

<?°(l,..,n | n+l,..,2n) = (-l)"(T/ 1 .../ ( J ti+1 ...7 2n ) is 
straightforward, i.e., the right-hand part of this relation 
will contain n! terms of ordinary Wick type (chain dia- 
grams) and also the different products of irreducible func- 
tions with the same total number of operators. These def- 
initions are the simplification of the corresponding ones 
for Hubbard and other lattice models proposed eaarlier 
by us. The calculation of the simplest irreducible func- 
tion like, for example, g^ lr (12|34) is rather cumbersome 
but straightforward. Then, we determine the Fourier rep- 
resentation of g^ lr [a i , Ti ; CT2 , T"2 1 o"3 , T3 ; o"4 , T4] , for which 
exists spin and frequency conservation, a\ + 02 = 
kd + u>2 — — UJ4 — 0, Thus, we dispose of the rules to 
deal with chronological averages in thermodynamic per- 
turbation theory. 

First of all, we determine the thermodynamic potential 
F of the system, 



F = F -^ln(U((3)) , 



(7) 



where Fq refers to the free impurity atom and conduction 
electrons. The diagrams which determine the thermody- 
namic potential have no external lines and represent the 
vacuum. Figure [T] shows the simplest connected vacuum 
diagrams of the normal (non-superconducting) state with 
the process of hybridization of C- and /-electrons. The 
zero order propagators of conduction and impurity elec- 
trons are represented by solid and dashed lines, respec- 
tively. These lines connect the crosses which depict the 
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FIG. 1: The simplest connected vacuum diagrams of the nor- 
mal state. The diagram (a) is of second and (b) and (c) 
are of fourth order in the hybridization interaction. Here, 

Vn = Vn/VN- 



impurity states. To each of the crosses two arrows are 
attached, one of which is ingoing and the other outgoing. 
They stand for annihilation and creation of electrons. 
The index n denotes (a n , r„) for impurity and (k n , <r„, r n ) 
for conduction electrons. The rectangles with 2n crosses 
represent the irreducible gn^ r Green's function, g^ lr . 

Besides the vacuum diagrams of fourth order in Fig. Q] 
(b) and (c), there is also one disconnected diagram com- 
posed of two diagrams of type (a) containing an addi- 
tional factor 1/2!. This is repeated in high order of per- 
turbation theory and permits us to formulate the linked 
cluster theorem, which has the form 



<£/(/3)) = exp (U(P))l 



(8) 



where (C/(/3))g contains only connected diagrams and is 
equal to zero when hybridization is swichted off. 

Now, we consider the diagrammatical analysis of the 
perturbation series for the renormalized propagator in 
Eq. (3). The simplest contributions to g in second order 
are represented in Fig. [21 All diagrams contain two exter- 
nal points with attached arrows fixed by the arguments 
of the Green's functions. At the inner points of the dia- 
grams a summation over <r n ,k n and integration over r„ 
is included. The new diagrammatical elements which ap- 
pear, are the irreducible two-particle Green's functions. 
The renormalization of these quantities is not considered 
here, only the renormalization of the propagators due to 
hybridization. 
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FIG. 2: The second order perturbation contribution to the 
impurity electron propagator. 

The final equations for the renormalized Green's func- 
tions are more conveniently written down in Fourier rep- 
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resentation. For the propagator of conduction electrons 
we obtain: 

G CTCT ,(k,k'M = <w<WG°(kM 

+ ^j^G* (kM W G>' (k'| iu;). 

(9) 

This renormalized propagator is expressed by the full 
propagator g of the impurity electrons. Now, it is neces- 
sary to obtain the corresponding equation for the full im- 
purity electron propagator. Because the subsystem of /- 
electrons is strongly correlated, we have to introduce the 
correlation function Z aa i which is represented by strongly 
connected diagrams containig irreducible Green's func- 
tions 0,0,11] and k a (iw) = g°(ioj) + Z a (iw). To this we 
introduce new quantities: 



E ^^(kr.^nlk^,^) = ^El^J 2 

kik 2 ki 

x G° lCT2 (k|n - r 2 ) = 5 aia2 G ai ( Tl - r 2 ), (10) 



which allow to simplify the structure of equation for the 
/-electron propagator, which can be written as 

g„(iu) = {A a (iuj) 

- Gl{-lUj)K a {lUj)k W {-iL0)}/d a {iL0). (11) 

d a (iu) = [1 - k a {iu)Gl{iu)][l - K w {-iu)Gl{-iuj)]. 

Equation (9) has originally been derived by Anderson [l[ 
by using the equation of motion method in which the 
propagator g a (iuj) corresponds to the t-matrix for non- 
spin-flip scattering. 

The main approximation which we employ here corre- 
sponds to the ladder approximation (which contain the 
main contributions of spin and charge fluctuations) and 
is schematically shown in Fig. [3] We investigate the prop- 

G(k 2 [T2T2 | kioyri) 



g 2 (fir; (Tin | ani\ a r ) 



FIG. 3: Diagrammatic representation of the approximative 
scheme used in the evaluation of the correlated function. The 
solid double lines with arrows depict the renormalized one- 
particle Green's function of conduction electrons. The rect- 
angle depicts the irreducible Green's function of impurity elec- 
trons. 

erties of the normal state of the single-site Anderson 
model by using Eqs. (9) and (11) for the corresponding 
Green's functions of conduction and impurity electrons. 
In the ladder approximation only the simplest irreducible 
Green's function g^ %T is iterated many times. In the 



paramagnetic phase (a' = a) the equation for Z reduces 
to 



Z a {iu}) = ~ E Gff i 



x g^' ir [a,iu;ai,iivi\cri,iuji;a,icj], (12) 



where 



G ^ IUJ ) = Jf E K^G^ksM. (13) 

kik 2 

Using Eqs. (9) and (11) G can be written as 



G a (iu>) = G° a (ico) + [G°M]VM 



(14) 



Then, using the definition of the correlation function, 
k a (iu), and approximation (12) for Z ff (iu>), we obtain 
the final integral equation for k a : 



a ct m = fl c?>M-±5:£i 



gf )lr [a, iu\a-L, iu 1 \ai , iu\ ; a, iw] 



(15) 



For the symmetric case when e/ = —U/2 < and 
ef + U = U/2 > 0, we have more simple equations: 

g°(iuj) = iu > {(iuj) 2 -(U/2) 2 }, 
Z = 2(1 + exp(PU/2)), n a = 1/2 (16) 

together with the antisymmetry property of the zero- 
order impurity Green's function, g^(— iu>) = —g„(iu>). 
In addition, if we assume that the matrix element V(e) 
and the bare density of states, po(e-) are even func- 
tion of its argument, then G„(iw) is also antisymmetric, 
G a a (— iuj) = — G a a (iu)). This allows us to search for the 
antisymmetric solution, 



k a (—iu) = -k a (iw), 



(17) 



of Eq. (15). The analytical continuation of G®(iw) has 
the form 

G° a (E + i5)=I(E)-iT(E), 

\V(e)\ 2 



I(E)=pJ<kp (ey- E _ t 
T(E)=7rp a (E)\V(E)\ 2 , 



(18) 



where 1(E) is the principal part of the integral. This 
function is antisymmetric. P(E) is the band width of 
the virtual level and an even function of the energy. The 
symmetric impurity Anderson model is easier to treat 
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since it leads to a simple form for the irreducible two- 
particle Green's functions in the different spin and fre- 
quency channels yielding 



A a (ico) 



3(iV2TOM 



+ 



[(zlo) 2 - (U/2) 2 ] 2 [l - A a (uu)G°(iLo)] 
iu) 



(ia J ) 2 -(U/2) 2 ' 



(19) 



where we have used results obtained for the one-band 
Hubbard model @, 0, HI ■ There are two solutions of Eq. 
(19). The physical one with the correct asymptotic be- 
havior, h a {iu>) — > 1/iui when \u>\ tends to infinity, has the 
form: 



A a (iuj) 



1 



2G°(t W )[M3-(E72)s 



{[(iu) 2 -(U/2) 2 + iLoG° a (iu)] 

[(iui) 2 - (U/2) 2 - iu}G°(iaj)] V1-12QM}, 

(20) 



where 



Q(iu) 



(U/2)G»(iw) 



{iuj) 2 -(U/2) 2 -iujGl{iw) 



(21) 



Finally, we obtain for the renormalized impurity electron 
propagator: 

g(iu) = 2^g a (iu)/G°(iuj) + 3Q(iuj)[g a (iLu) 

-1/G°M] 2 }{[1/G°M 

-g°(iuj)}[l + ^l-12Q{iLu) 

+ 6Q(iu}){g (iu})G°(iu)-l}}\ . (22) 

From the analytical continuation we obtain the spectral 
function of impurity electrons, which is given in Hubbard 
I approximation by 



A){E) 



2T(E) 



[E - (U/2) 2 /E - 1(E)] 2 + T 2 (E) ' 



(23) 



Equation (23) has two resonances at E = ±Eq, where 
Eq is determined from the condition that the brachet in 
the denominotor of (23) is equal to zero. These reso- 
nances are broadened by T(E) which determines height 
and width of both resonances. This approximation yields 
a vanishing spectral function on the Fermi surface, where 
E = 0. Only beyond Hubbard I approximation, when 
strong correlations are taken into account, we obtain 
from Eq. (20) the spectral function of impurity electrons 
in the vicinity of the Fermi level having the following 
form, 



ME) 



i2(2r/c/) 2 /r 



1 + 6(2T/U) 2 + yft + 12(2r/[/) 2 ' 



(24) 
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FIG. 4: Spectral function Af(E) x W for different values 
of the theory parameters as function of energy E/W in the 
Hubbard I approximation (case (a)) and in the our ladder 
approximation (cases (b) and (c)). In the cases (a) and (b) 
U = leV, W = 2eV, F = 0.1 eV and in case case (c) 
U = 1.5eV, W = 2eV, r = 0.15eV\ 



In conclusion we state that our diagrammitc approach 
allows to describe in ladder approximation for the ir- 
reducible Green's function two resonances at energies 
E = ±Eq and the feature at E = given by Eq. (24). 
The distance between two peaks (Figure @| is determined 
by the value of U but the height and the width of the 
peaks by T. Beyond the Hubbard I approximation the 
non zero values of spectral function appear for E = 0. 
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